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ABSTRACT 


A theoretical method of predicting the performance of subsonic, axial 
flow, multistage turbines is presented together with the digital computor 
program for computing all the dimensionless performance parameters required 
to completely define turbine performance. A small two stage turbine, for 
use in space vehicles, was used to demonstrate the application of the method. 
A complete set of performance maps were drawn and analyzed. The dimensionless 
performance parameters for any given flow condition could be obtained from 
the maps. 

The computor program proved to be extremely flexible and useful. The 
effect of blade row redesign could be easily determined. Comparison of the 
extremely limited amount of open cycle test data with program results showed 
that the method would provide a design engineer the means of predicting the 
performance of a given turbine design. The accuracy of such a prediction was 
shown to depend greatly upon the estimation of rotor tip clearances and the 
measurement of flow areas corresponding to the clearance. The computor pro- 
gram provides a means for trial and error determination of the rotor tip 
clearances when operating at high temperatures if accurate test data is avail- 
able. 

Although the computor program was written in Fortran language for the 
Control Data Corporation 1604 Computor at the U.S. Naval Postgraduate School, 
it should be compatible with most computor installations through out the 


Country. 
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SYMBOLS AND UNITS 


DEF INITION UNITS 
Gas flow area measured normal to flow direction So Lit. 
Minimum blade passage gas flow area sq. in. 
Blade opening or distance between blades at the in. 


minimum cross section 

Lift coefficient based upon vector mean velocity 

Theoretical velocity for isentropic expansion from ft/sec 
stagnation pressure at entrance to static pressure 

at discharge 


Specific heat at constant pressure Btu/1lb a 


Coefficient of drag on blades created by tip 
clearance pressure losses 


Coefficient of drag on blades created by secondary 
flow pressure losses 


Blade chord in. 
Diameter ine 
Gravitational constant B27 Gs epeeer 
Total absolute enthalpy Btu/1b 
Total relative enthalpy Btu/lb 


Shape parameter; Energy thickness 


Total enthalpy drop Btu/1b 
Static enthalpy Btu/1b 
Blade height in. 
Inside diameter of turbine annulus msi. 


Incidence angle of flow onto a blade row, given 
by the difference between gas flow angle relative 
to blade inlet and blade inlet angle 


Stalling incidence 


Conversion factor Jisatt—ib/Btu 


Vaid 





Radial rotor tip clearance 

Rotational speed 

Polytropic exponent 

Outside diameter of turbine annulus 

Total pressure 

Total pressure drop 

Static pressure 

Free stream dynamic pressure 

Gas constant 

Entropy 

Blade spacing 

Gas total temperature 

Equivalent total temperature at rotor inlet 
Equivalent total temperature at stator inlet 
Gas total temperature drop 
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Maximum blade thickness 
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Relative velocity 

Tangential component of relative velocity 


Gas mass flow rate 
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yee 


rpm 


Tie 
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psia. 

ib/ét? 
ft 1b/1b °R 
Btu/1b °R 
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il 
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ft/sec 
ft/sec 
ft/sec 
ft/sec 
ft/sec 


lb/sec 
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Stator gas flow angles Degrees 


Stator blade angles Degrees 
Rotor gas flow angles Degrees 
Rotor blade angles Degrees 


Ratio of specific heats 

Small finite interval 

Shape parameter; Displacement thickness 
Expansion loss coefficient 

Rotor tip clearance loss coefficient 
Profile loss coefficient 

Secondary loss coefficient 

Total loss coefficient 

Trailing edge loss coefficient 
Efficiency 

Flow area reduction factor 


Power coefficient; also a factor defining secondary 
losses 


Gas density sited ee 


Flow function 


Ww) Angular velocity rad./sec 
SUBSCRIPTS 

er Crieieat 

D Diffusor 

e exit section 

ed. equivalent 
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N Nozzle 
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5 stator 
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e Total 

A Trailing edge 

u Tangential 
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Note. See Appendix III for Fortran names and their meaning. 
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A METHOD OF PREDICTING THE PERFORMANCE OF AXIAL FLOW TURBINES 
USING A DIGITAL COMPUTOR TO DEVELOP PERFORMANCE MAPS 


INTRODUCTION 


Turbines for driving auxiliary equipment on rocket and space vehicles 
have become on of the most exacting pieces of turbomachinery in use today. 
The premium placed on space and weight in these vehicles requires that the 
size of a turbine be small, yet it is necessary to extract the maximum work 
for a given input. This means that the turbine must be designed to have 
high specific work output and efficiency. The effect of size on the effici- 
ency of turbines has not been fully determined, however small turbines re- 
quire that clearances be small, flow areas exact, leakage losses low, and 
deflections of the flow in the rotor blade less than 105°, Ref. 1. Regard- 
less of size, calculation of the performance of a turbine must be made using 
approximate methods. 

Although an exact analysis of viscous, compressible flow through an 
axial turbomachine can never be made, the demand upon the design engineer 
for an accurate estimation of the performance of a given design has become 
greater. The greatest difficulty in making an accurate prediction is con- 
tending with the large number of variables that play a role in the overall 
performance of a turbine stage. Since most engineering firms have access to 
electronic digital computors, a method of analyis involving their use was 
suggested. 

Manual calculation of the efficiency and work output of a multistage 
turbine for one set of operating conditions is tedious and time consuming. 
For a two stage turbine, approximately six man-hours of work are required 
to obtain one constant speed characteristic point as defined by the flow rate, 
overall pressure ratio, efficiency, and power coefficient. Many combinations 
of the input parameters are necessary in order to make an accurate estimation 
of overall turbine performance. In order to provide a quick and accurate 
method of evaluating the performance of subsonic, axial flow, single or 
multistage turbines, equations were developed and programed for computor 


solution so that performance maps could be drawn. 





Under the assumption that a multistage axial flow turbine design pre- 
sented for evaluation would give the blade shapes, flow angles, diameters, 
clearances, spacing, and all other pertinent data, the computor program de- 
veloped will allow an accurate evaluation of the performance to be made. 

The program has been restricted to subsonic flow conditions. No extension of 
this program can be made to include turbines operating in the supersonic 
range since a complicated iteration process would be necessary in order to 
determine the pressure ratio across each blade row. To preserve flow con- 
tinuity, the deflection of the outflow from a blade row would also have to 


be considered, Ret. Ze 


II. Method in General 
A. Simplifying Assumptions 

Once a turbine design has been formulated, any estimation of the 
performance must be made using three dimensional flow conditions as a basis, 
Since three dimensional effects are so important. Due to the complexity of 
such flows certain assumptions have been made for simplification. 

The flow has been assumed to be adiabatic, steady, turbulent, and 
axial at entry to the nozzle blade row of the first stage. Frictional forces 
have been ignored in the region between the blade rows, since velocity gradi- 
ents in that region are much smaller than those found in the boundary layers 
along the surfaces. The flow between the blade rows has been considered to 
be axisymmetric and steady, depending solely upon the conditions imposed by 
the blade rows ahead of a given region. Interference effects between the 
rows of blades have been ignored. At the mean radius the flow was considered 
axisymmetric, having axial and tangential velocity components. Since the 
annulus radial height is usually small compared with the mean radius, the 
changes of the flow in radial direction have been ignored and the mass flow 
rate was determined for the conditions at the mean radius. 

An assumption was made that the gas outflow angle from a blade row 
is independent of incidence and Mach number. While this assumption is not 
precisely true, no appreciable error will be incurred over the efficient 


operating range of a subsonic turbine, Ref. l. 





Ot course the evaluation of the performance of any turbomachine, 
no matter what method is employed or how many place accuracy a computor can 
achieve, depends primarily upon the accuracy of the estimation of the losses. 
These losses are imposed by the frictional forces in the boundary layers along 
the blade surfaces, by mixing, and by clearance effects. The major simpli- 
fication of considering the flow path through each stage at one diameter only 
requires the further assumption that in any one cross-sectional plane of the 
flow between adjacent blade rows the total pressure, total temperature, and 
axial velocity are the same at all points. Such an assumption, though widely 
divorced from fact, may yield correct overall characteristics of a stage if 
the loss coefficients used are equal to the momentum mean values over the en- 
tire cross-sectional plane. 

Since the accuracy of the performance calculations rests mainly 
upon the loss coefficients, it is imparitive that they be as accurate as 
possible. American and British design methods and cascade test data were re- 
searched in order to determine a basis for loss determination. American 
methods of determining blade row losses are based mainly upon theoretical con- 
siderations. A British method of predicting the loss coefficients, Ref. 3, 
based upon test data derived from overall tests on a variety of turbine 
Stages having blading approximately midway between impulse and reaction types 
was chosen. The experimental data from which the loss coefficients were ob= 
tained had been deduced from tests made at Reynolds numbers in the range of 
x 10? to 3 x Te. therefore the method could only be used for turbines 
Operating in that range. The test data apply mainly to blades having a con- 
ventional profile shape. Most blade shapes in current use in gas turbines 
fall within that catagory. The method allows the loss coefficients and 
effux angles in any blade row to vary with gas flow conditions and the angle 
of incidence. The loss coefficients are assumed to be uninfluenced by Mach 
number. This assumption is unlikely to cause an appreciable error unless 
the blades involved have a high degree of curvature on the upper surface near 


the trailing edge. 


B. Use of Dimensionless Parameters 


In order to simplify the analysis and to present the complete per- 





formance of any turbine graphically, the performance can best be determined 
using dimensionless parameters involving the different variables which in~ 
fluence turbine behavior. In this way the complete performance of a turbine 
under a variety of inlet conditons and speeds can be presented on three dia- 
grams. The four main dimensionless parameters used were refered flow rate, 
refered rpm, overall pressure ratio, and a power coefficient. The conditions 
at any point in the turbine have been refered to inlet conditions. Ref. 4 
presents a complete development of these parameters using Riabouchinski's 


theorem. 


1lII. Description of the Turbine to be Analyzed 

The turbine chosen for demonstration of the method and program was a 
small two stage axial flow turbine which was being evaluated for use in a 
Space vehicle. It had proven unsatisfactory in the few preliminary tests 
conducted, due to low efficiency and power output. It was believed that 
poor agreement between test results and calculated flow quantities was due 
to inaccurate calculation of the required minimum flow areas between adja- 
cent blades, and unrealistic assumptions of the loss coefficients. Small 
differences in blade thicknesses and blade angles between design drawings 
and the manufactured product can make large differences in the minimum flow 
areas. 

Fig. 1 shows a 5:1 scale drawing of the meridinal blade passage of the 
turbine. Included on the drawing is the location of the local minimum flow 
areas and cross sections of the blades. Fig. 2 illustrates the system adop- 
ted for defining the geometry of a blade row and the gas angles relative to 
a blade row. Table I gives the dimensions and angles of the blades and other 
pertinent blade data. All blade data were obtained from Fig. 1 and the blade 
section drawings, Fig. 3 through 9. Not all of the blade section drawings 
are presented. Section drawings of the Stator and Rotor II are included to 
show the actual blade shapes, the minimum distance between blades, and to 
allow determination of the angles and dimensions at the mean blade height. 
Values for a, t, ee s, Q@*, B*¥, c, and A, were obtained from the drawings 
for the mean diameter of each row of blades. Fig. 10 shows the blade angles 


for each row of blades and the turbine station designations. 








The tlow angle at the discharge of a blade row was determined from the 


empirical relation 


=] a 
Ot : = ¢o& 
exit S - (t/cos ax 


= 
This equation is widely used in Europe and the U.S.S.R. and is believed to 


be more accurate than the relation 


commonly used in the United States. The sign convengtion used designates 

all gas angles positive if the tangential component of velocity is in the 
direction of rotor motion. It should be noted that the angle of incidence 

of the flow into the nozzle blade row is 1S , Fig. 3. A mMoedtficatien of 
the nozzle blade design is contemplated which would reduce the angle of in- 
cidence to zero and reduce the losses. The trailing edges of the blades are 
relatively thick in order to prevent burn off at high temperature operation. 
All the blades converge to a minimum area shown on the blade section drawings. 
The area was determined by multiplying the average height of the minimum flow 
area between blades by the height of the blades. 

Nitrogen gas was specified as the working fluid. The specific heat 
ratio for nitrogen is a function of temperature and varies from 1.3/2 to 
1.347 over the temperature range of 780 to 1260°F, Ref. 5, however the 
variation was considered small enough to allow an average value of 1.36 to 
be used in all calculations. The average value of the gas constant over the 
temperature range is 55.16 £t/OR. 

The design point of the turbine was reportedly 18,000 rpm and 1200 ie 
In order to check this point, the operating characteristics of the turbine 
were investigated over the range of rpm between 10,000 and 19,000 at teme 


peratures between 780 and 1260 ye 


IV. Development of Flow Function Formulas 
In order to determine turbine performance for a given set of inlet and 


speed conditons, the course of an element of mass of gas was followed from 


one blade row to the next, with calculations performed at each station. 
In this way the basic performance parameters were obtained for a given 
refered flow rate and refered rpm. 

For steady adiabatic flow conditions the stagnation enthalpy along 
any streamline remains constant, both for absolute and relative flows. 
The mass flow rate for a given set of inlet conditions is constant and 


can be expressed by 


O 


Ww = 
\ RT | 


The polytropic exponent was used in the power terms because the flow 
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was adiabatic but not frictionless. See Appendix I. 

Friction within a row of blades reduces the overall efficiency. In 
determining the flow rates with friction, the rate of flow is governed 
by the area of the minimum cross section. The flow area is reduced by 
the build up of the boundary layer along the profile from the leading 
edge to the point on the blade corresponding to the minimum cross section. 
The amount of boundary layer growth depends primarily upon the blade pro- 
file and the Reynolds number of the flow. In this paper the result of 
this reduction of flow has been termed a loss and designated ( |. Dies 
loss coefficient influences the flow rate and is believed to have more 
effect on the overall performance of a blade row than the total loss 
coefficient, ( The efficiency of a blade row can be expressed in terms 


Ss 


of the temperature ratio dT/dT, om 
is e 
= qT/dT, = 1 te (2) 


The polytropic exponent can also be defined in terms ae 


n= B/(C (¥- 1) +0) (3) 


The development of this expression is shown in Appendix I. 





The mass flow rate can be expressed in the form of a nondimensionless 


flow function by rearrangement of Equation (1) 


| n+ 1 =<; 
PP) n | (3) 








¥ thy 2/n ,1 
2 on Pop) 


Since n is a function of the expansion loss incurred from entrance to the 
minimum area of the blade passage, ¢ is a function of Q , on and the 
pressure ratio P /P: For a given value of %, a plot of $ versus P/P 
will produce a curve as shown in Fig. 11 for each value of a The data 
for a complete set of such curves for values of C | fromweo LO .25 aden 
from 1.20 to 1.40 was obtained by Vavra using the CDC 1604 computor. The 
data for a specific heat ratio of 1.36 is given in Table IV. Using this 
data it is possible to determine values of P /P for calculated values of 
$ by two way interpolation. The table is useful in manual calculations. 
In order to adopt the flow function formula to a computor solution 
of the pressure ratio that would be accurate and minimize the time of 
computor utilization, it was necessary to calculate an approximate press- 
ure ratio by expressing the terms in binomial series form. Substitution 
of the first two terms of the expansion in the equation for > was made in 
order to obtain Ap/P . as a funce) omeol on sane > As a first approxi- 
mation P/P = Ee Gee - Ap) = 1/(1 - p/P) or p/P. = l ~Ap/P.. The 


complete development of the equation for the approximation 


a nr ee 

(P /p) = 1/(1L - n/3 CL Yh - 30% - LYSE - eee )) (6) 
is given in Appendix 1. 

Using this approximate pressure ratio, a value of ¢$ could be calculated 
and compared with the known value of b . By increasing or decreasing the 
approximate pressure ratio by an increment as necessary, the value of P /P 
that corresponded to the known $ could be obtained. This method of finding 
the pressure ratio corresponding to a given ? was made into a subroutine 


called "Ratio" for the computor program. 





V. Methods used in determining Loss Coefficients 

In order to calculate the changes in the pressure and temperature from 
point to point saeich the turbine, it was necessary to establish the pres- 
sure losses that were involved. The overall pressure loss occuring in a 
blade row was subdivided into a number of component losses which are de- 
pendent upon various variables that define the aerodynamic form of the 
gas flow and the geometric form of the blade row. These losses are de- 
pendent upon the angle of incidence of the flow into the blade rows. The 
component losses considered were 

a) Profile losses-~losses due to skin friction which causes the build 
up of a boundary layer on the blade profile. 

b) Secondary losses--losses resulting from non-uniformity of the three 
dimensional flow through a row of blades mainly caused by the interaction 
between the blade ends and the boundary layer on the annulus walls. 

c) Tip clearance losses--losses due to leakage of gas over the ends of 
the blade tips. 

d) Trailing edge mixing loss--a loss caused by the thickness of the 
trailing edge of the blades. 

The profile losses were determined using the methods presented in 
Ref. 3. The profile loss for a given blade row is first determined for 
inflow at zero incidence. The stalling incidence of the blade row was then 
determined, stalling incidence being defined as the incidence at which pro- 
file loss is twice the loss at zero incidence. Profile losses at incidence 
other than zero were based upon the assumption that the ratio of profile 
loss at any incidence to profile loss at zero incidence is a function of the 
ratio of incidence to stalling incidence. 

The profile loss coefficient at zero incidence was assumed to be a 
function of the discharge flow angle, the ratio of the blade angle to the 
discharge flow angle, the pitch to chord ratio, and the thickness to chord 
ratio. Ainley, Ref. 2, defines this loss coefficient as Y = loss of total 
pressure divided by the total pressure at discharge minus the static pres- 


sure at discharge. 


Se re 





For this equation to be valid the total pressure at exit would have to be 
measured far downstream of the blade row. The pressure would also have to 
be an average value. The total pressure at discharge in the immediate 
vicinity of the trailing edge of the blade is known to fluctuate. It was 
considered more correct to express this loss coefficient in terms of Ah, | 
rather than in terms of OP/(0/,V>). 

Using the differential form of the Energy Equation, where dq = 0 for 


isentropic flow, the loss of a blade row can be expressed in terms of Ah 
dq = du # p dv = dh - v dp = 0 


dh = v dp = RT dp/p = dp/f 
or Ah Ap/? 


lI 


{ 


Substitution of the expression @Ah for Ap in Equation (7) gives 
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as shown by Fig. 12. The right and left sides of the equation are not 
exactly equal since the temperature at inlet to a blade row is not the 
same as at the exit. However this difference is so small for an individual 


blade row that it can be considered insignificant. 





y= 
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Expansion Process for a Blade Row 





By rearrangement 











2 v2 
See _ DS 
ne 2e5 Cite toss OMMEOSS = Toy “teal (8) | 
and aT can be considered a new loss coefficient, ( ; 


The imperical equations for determination of aoe Yo» and Yy Ait. 6 

incidence were obtained from Ref. 6 and are presented in Appendix []I. Com- 

Plete calculation for the loss coefficients of the Nozzle and Rotor I are 

also presented in Appendix II. The calculations for the loss coefficients 

for the Stator and Rotor [II are similiar to those for the Nozzle and Rotor I 
respectively. Variations of the losses with gas inlet angle to a blade row 

at large positive and negative incidences are uncertain, but according to | 
Ainley, Ref. 6, reasonable correlation of test and calculated turbine per- 
formance has been obtained by restrictingthe use of the equation for M. and 
YE to the range of -1.5 <i/ii¢ 1.0. At values of i/i_>1.0 the secondary 
and clearance loss coefficients should be assumed constant and equal to the 
value when i/i. = 1.0. Similiarly when i/i, <-1.5 the values of X. and Y, 
equal to those obtained for i/i. = ~1.5 should be used. 

The profile loss coefficient a for conventional section blades at 
zero incidence can be obtained from Fig. 13. Values of positive stalling 
incidence of cascades are shown in Fig. 14. The variation of profile loss 
with incidence is given in Fig. 15. Using the curves presented in these 
Figures, all of which were reproduced from Ref. 3 and 6, values of C at 
ten degree intervals of incidence angle were determined and are presented 
for each blade row in Tables V, VI, and VII in Appendix II. 


In determining the secondary loss coefficients Ainley made the basic 


assumption that 
2 
= S 
c= Ac," /(8/c) 


where ee primarily dependent upon the degree of acceleration imparted to 
the gas as it flows through a blade row. A similiar assumption was made in 
regard to the tip clearance loss coefficients in stating that the drag coe 


efficient can be expressed by 


C., = Constant (k/h) ¢, “/ (s/c) 


Dk 


10 





The value ere used in determining Ye was obtained from Fig. 16. It should 
be noted that secondary and tip clearance losses in a blade row having fixed 
inlet and outlet gas angles are independent of s/c, thus optimum pitching 
of a row of blades is obtained using the pitch that gives the minimum pro- 
file loss. 

The equations presented in Ref. 6 for calculation of the trailing edge 
loss were purely theoretical and were not considered to be as accurate as 


the equation 


bce ee Cp (9) 


) Ox/a 


developed by Markov, Ref. /, using average measured values of the shape 
Parameters. Vavra, in an as yet unpublished paper, Ref. 8, explains the 
development of equations for ( - and OF The assumption that the 
boundary layer thickness is Ge same at the throat as at the trailing edge 


was made, Fig. I1/. 





Fig. 17 


Boundary Layer Thickness Diagram 


fe 








Such an assumption would mean that further expansion of the boundary layer 
between the minimum cross sectional area and the trailing edge does not take : 
Place. Soul has been defined as the sum of the displacement thicknesses of 
the boundary layers on both sides of a profile at the trailing edge of the 


blades. The flow rate was expressed as 


ies 


W = PAV, =PyVo(a +> = 5V oC) = eee 


or w= Poona Y 


where ]?/ is a factor which accounts for the reduction in the flow area at 

the minimum cross-section of the flow passage between blades. Vavra developed 
; > e 

a theoretical expression for, from curves of - >—, versus the pressure 


1 - 


ratio across a blade row. 
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ao ~ f Sarat Sie 

1 -) discharge 
For a X of 1.36 

S 
/s 
C. =.9 1-95") (10) 


Energy thickness was defined as 


ane aes 
H Oa 
kik 
and a value of H = 2.2 was given in Ref. 7. (1 -C)) was expressed as 


a ADs ii 
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Upon rearranging terms 
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from which the following equations were obtained 


Z 


Ue TET, Oey 
Oe a (12) 


These equations were utilized in calculation of the trailing edge loss and 
the reduction in flow area due to the build up of the boundary layer along 
the profile from the leading edge to the throat of the blades. 

The values of the total loss coefficients given in Tables IV through 
VI are for the design tip clearance of .033 in. for Rotor I and .021 in. 
for Rotor II. These clearances are reduced when operating at elevated 
temperatures. Since it is impossible to say exactly what the clearances 
are at the operating temperatures, loss curves were also drawn for tip 
clearances of .005, .010, .015, and .020 in. The loss curves for each 


blade row are presented graphically in Fig. 18, 19, and 20. 


VI. Development of Basic Equations 

The equations used in calculating the state of the flow at a given 
point in the passage through the turbine were based upon the equations 
for absolute and relative flows. These equations have been fully developed 
in Ref. 1 from the basic Equation of Motion, Equation of Continuity, and 
the steady flow Energy Equation. For steady adiabatic absolute flows 
through stationary blade rows the total enthalpy along a particular stream- 
line is constant. Likewise for steady adiabatic relative flows through 
moving blade rows the total relative enthalpy along a relative streamline 


is constant. Expressed in terms of velocities and static enthalpies 


v2 
H=h+>— + ez (13) 
mer absolute flows and 
2 22 
W WwW R 
Hp = h + 7 9 + gz (14) 


for relative flows. The effect of the gz term has been considered negli- 


gible. Ref. 1 also explains in detail how Euler's Turbine Equation 
| Ueaey ae 
ul 2 ue 


1 
An = C, At, = a5 
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can be developed either by combining the enthalpy equations for absolute 
and relative flows between two points or from the general Momentum Equation. 
While this relation has not been included in the computor program, it can 
be used as a check on the solution of the total enthalpy drop through the 
turbine if manual calculations are attempted. 

The equations used in the computor program were developed and ar- 
ranged so that calculations would be repetitive for any number of stages. 
The nozzle was considered to be the same as any other stationary blade row 
except that the velocity vector of the inlet flow was considered to be axial 
in direction. For each blade row the absolute velocity, relative velocity 
in the case of rotor blade rows, was first determined using the losses ob- 
tained from a knowledge of the inlet flow. The temperature at discharge 
from a blade row could be calculated from this velocity. Also the velocity 
triangle can be drawn if the velocity and angle of the inflow are known. 
For each blade row the flow function method was used to find the pressure 
ratio across it. All losses were considered in percent of the theoretical 
kinetic energy at discharge from a cascade. All pressures at inlet toa 
blade row were expressed as a ratio of the total pressure at that point to 
the pressure at inlet to the turbine. 

A form of the basic equation for steady adiabatic flow was formulated 
whereby an equivalent temperature and pressure at the minimum area of the 


blade passage was determined. The flow function could then be expressed as 





 - y IT, Teq/to R (16) 
i. P P /P g A 
O eq oO Os 


where A, was the minimum flow area for the cascade. 
A. Nozzle and Stator 
The equations for calculation of the thermodynamic flow process 
Occuring in the Nozzle or Stator are based upon the fact that the total 
enthalpy along a streamline is constant. For steady adiabatic flow, the 
enthalpy drop at constant entropy can be expressed as 
2 


V 


: lis 
On, “sj — ah 





rae (loss coefficient x Ah, .) 
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or in terms ot velocities 





v2) Seemed ie 
Lis) tae ne lis (17) 
2gJ wee N 2gJ 


Sance h = T, the relationship can also be expressed as 


ger 7 eee Svar, cre) 


and shown graphically in Fig. 21. 
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Fig. 21 


Temperature Change for a Blade Row 


The velocity at the discharge of a stationary row of blades was 


determined from the following relations 


v,* ( 
——— = = ~ 4 
2gJ au me N ye 
= C, 4T = C, (AT. - y OT; ) 
-¢, 47). Ly (19) 
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AT, . can be expressed in terms of pressure ratio 


AT. =|1 eae 
1S 


When a substitution for AT,. is made 


v." p oe 
pS earre t,t; [bg > | 


For simplicity and ease of calculations all velocities were divided by 


FOO: <s6 


, Sel 


V 
= 5.007 (-G.) [4 “Gh Fy bea so .007 Bi ot [Gg soo }20) 


tj 
100 








When the pressure ratio P,/P)> corresponding to the value of the flow 
function $ is substituted into Equation (20) the value of Vy can be 
calculated. The equation is in the form of 

Z 


vi ( 
= constant (1 - nN AT, . 


100 








Bince I -T, = (1) - Cn 


2 
“1 (21) 
100 = constant (T “ T)) 





Equation (21) was used to determine the temperature at station l. 


B. Rotor 
All values of the velocity triangle representing flow conditions 
at inlet to a rotor blade row, Fig. 22, can be calculated when the peripherial 
speed of the rotor blades, the absolute velocity of the inlet flow, and the 


angle of incidence are known. 
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Velocity Triangle, Rotor Inlet 


The following equations allow ease of calculation of the magnitude and 


direction of the vectors: 


vol = vy sin a, 
vq = va cos a 
W = Vv - U 


-1 
2 eae ine ae 
If the loss curves for the first rotor are entered with the inlet flow 
> - . ’ 
angle By> the losses (, and © 4 can be obtained. In order to find the 
static pressure after the rotor the flow function must be based upon relative 
flow relations. In the rotating blade rows the relative total enthalpy must 


remain constant along a given relative streamline 


Wo ue 
Hp h + Ded Qed constant C22) 
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Since the radius of the mean flow path at rotor entrance differs from that 
at rotor exit, there will be significant difference in the peripheral speed 
at these two points. The peripheral speed at a point is given by the ree 
lation 


U = xND/720 


where the diameter is measured in inches. 
The total enthalpy at inlet to the rotor blade row is the same as the 
total enthalpy at the minimum area, therefore 


2 2 x 2 
h "l - mcs ee e ie a e Ekg 
eon 223 = 2 22] 2eJ 





* w 
The sum of hy and Wo /2gj3 can be considered to be an equivalent enthalpy, 
ec » Similiar to the sum of h and Tet for a stationary row of blades. 


When the variation in the peripheral speed is taken into consideration 


H = hy # (w # U ie U,*)/ 2gJ (23) 


eq. 2 


The equivalent enthalpy can aliso be expressed as oe Ae . An equation for 
determination of the equivalent temperature can be derived by dividing 


Equation 23 by 8 


2 
it = T, f (W, # U 


eq. 


2 2 
» ~ Uy) 2guC, (24) 


The equivalent total temperature and pressure at inlet to the rotor 


was designated TR and P. respectively. In ratio form 


Gout 
Pay/Py = (Tpy/T)) Bol Soa 


The flow function for the rotor can be written as 


— NTa \ fag. 
p R Pl eC a 


where A, is the area at the minimum flow cross section. By refering all 


temperatures and pressures to turbine inlet conditions, a nondimensional 


—- 
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equation is obtained for the \- function 


W a Vai od 


p/P, fo 


The relative velocity after the rotor can be expressed as 


W W + U 3 - U é 
ae |, tee eee ee 
Ps 1” 2eJ 2 
W ee; au 
or 2 -CT iL Zz 1 -cT 
283 jee 285 2 
= C coe - T,) (26) 
me. 1 -/ A (27 
ae” 2gJ GR Tis ) 


These relations are presented graphically in Fig. 23. 
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Temperature Change for a Rotor Blade Row 
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The velocity triangle representing flow conditions at exit from the 
rotor blade row can be determined from the peripheral speed, the relative 


velocity, and the angle of exit, Fig. 24. 





Velgelty lomanete, Rotor Exit 


The equations involved are 


Wi2 = W, sin Bo 


= W, cos B. 
Vu2 = U5 : we 
2 Z 2 
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The pressure ratio Poo/P at entrance to the next stator blade row can be 


calculated 
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The procedure outlined above can be repeated for as many stages as 
necessary depending upon the turbine design. Fig. 25 is a complete T - § 
diagram representative of the two stage turbine considered. Sample calcu- 


lations are presented in Appendix IV. 


C. Diffusor 
For stationary gas turbine power plants an efficiency for the 
Diffusor of 70% is commonly accepted for a flow that departs axially from 
the last row of blades. For a flow that is discharged from the last blade 
row at an angle to the axial direction, the actual diffusor efficiency was 


considered to be 


Ft ate (cos a)? 


The overall turbine efficiency was defined in terms of enthalpy 
nad AM 
roe. LPA 


1S S 





Defining turbine efficiency in this way accounts for recovery factors. 

The ability of the diffusor to transform the kinetic energy of the 
flow at exit from the last blade row is a function of the efficiency of 
the diffusor and the difference between the kinetic energy of the flow at 
inlet to and discharge from the diffusor. In terms of temperature change, 


this relationship can be expressed as 


~ erecoveny = Ni 





where the velocity of discharge was calculated using the Continuity Equation 


wt D Ay Vy = constant 
‘= W RT), 
ae Ty, 


By summing the work output of the individual stages the overall work 
of the turbine was obtained. The specific work output of the turbine is 


equal to the enthalpy change AH between the inlet to the first blade row 
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and the discharge from the last blade row. For the two stage turbine con- 


Sidered 


AH = (U, Ve ee V 2/85 + (U, ve a Vg) /85 = C40, 


The power output was desired in coefficient form. Since the power 
can be expressed as the mass flow rate times the enthalpy change, a suitable 


coefficient form is 


P w\T AT 
~ “- 1.055 
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KW 
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A= ste = ¢ (27) 
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where / T= (CH. - 4r/C,, Solent oe 
The overall pressure ratio of the turbine was determined by first 
calculating the ratio of the total pressure at discharge from the diffusor 


to the static pressure at entrance to the diffusor. 


os Dy NT y! 
—— 4 LSD Sa 
P./p, = (T,/T,) f-1 (aes x-1 (28) 


Using the pressure ratio across the turbine blade rows, P /Py? the overall 


pressure ratio was calculated 


ee Z (22) 


The isentropic temperature drop through the turbine could be found 


using the relation 
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Turbine efficiency can be computed using the definition 
n ; AH ; AT. 
io (AT 
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In order to compute the overall velocity ratio of the turbine the 
mean average diameter of the flow passage of the turbine was calculated. 
Since the theoretical velocity oe for isentropic expansion from a stag- 
nation pressure at turbine entrance to the static pressure at turbine dis- 


charge could be expressed as 


Cc. = \2g¢, AT 


an equation for the velocity ratio was 


U x N D 
avg. | avg, 


; = 
3 720 \28JC, yNerer, 


VIL. Computor Program 
A. General 

The Control Data Corporation 1604 digitial computor at the U.S. 
Naval Postgraduate School was utilized to provide rapid and accurate solu- 
tions to the turbine performance equations. The source program was written 
in the most basic and familiar version of Fortran language so that the pro- 
gram would be compatible with other models and makes of computors. In order 
to clearly document and visually present the step by step procedures of the 
program, flow charts were drawn for the main program and the major sub- 
routines. The flow charts and selected versions of the basic program are 
Presented in Appendix III. Transfers of control and test routines are 
shown more clearly by flow charts than if described in words. A table of 
Fortran names, equivalent symbols, and meanings is also presented in Appendix 
iL... 

B. Main Program 

The complete Fortran program was a composite of a main program ‘and 
several sub-programs. The main program was used for control and input-output, 
while the subprograms performed the repetative calculations. Input data such 
as turbine blade row dimensions, blade angles, and loss coefficients were 
placed in one dimensional arrays. All constants and variables which were re- 
quired in the main program and one or more sub-programs were included in a 
Common statement so that communication between the main program and sub- 
programs was possible. The values of diffusor cross sectional area, average 


mean flow diameter, specific heat ratio and gas constant for the working 
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fluid were considered constant and initialized in the main program for use 
in all sub-programs. 

Only the parameters necessary to demonstrate the performance of the 
turbine and allow development of the performance maps were normally printed 
out for each combination of refered rpm and refered flow rate. The para- 
meters most representative of turbine performance were refered rpm, refered 
flow rate, efficiency, power coefficient, overall pressure ratio, and speed 
ratio. Of course the print out of the solution to every equation was pos- 
Sible. Such a print out was made for the test case, Appendix IV. 

The performance of the two stage turbine considered in this paper was 
desired over the temperature range from 1240 to 1720 °R and the speed range 
from 10,000 to 19,000 rpm. The resulting range of refered rpm was 240 to 
540. A "Do" loop was inserted which allowed calculations to be made at 
any desired interval over the range. For ease of plotting and completeness 
of coverage, an interval of 50 was initially chosen and later reduced to 10. 
The performance maps show only the curves for refered rpms of 240, 290, 340, 
390, 440, 490, and 540 in order to allow curve separation and prevent con- 
fusion, however the smaller interval was necessary in order to accurately 
complete the performance maps. In order to compare the theoretical com- 
putations with actual test results a single refered rpm was programed so 
that the Do loop would start and stop on the same value. The computor time 
for calculation of one test run was approximately one minute and ten seconds, 
compared to two minutes and fifty eight seconds for complete coverage of the 
refered rpm range using an interval of 10. 

The performance of the turbine was desired at all values of refered 
flow rate between zero and that value which would cause a turbine blade row 
to choke. Although the flow rate Do loop was programed to start at .1 and 
continue to 5.0 in steps of .l or .0O1, the inflow angles to the blade rows 
were so great at the lower values of refered flow rate that the range of 
angles over which the loss curves were valid was exceeded until a flow rate 
of approximately 2.5 was reached. The range of the Do loop was reduced by 
Starting at a valwof 2,0 in order to reduce the computor time involved. 

A flag, ICR or IBR, was set to test whether the angles exceeded # or - 70° 
If the flow angle was excessive at inlet to either the rotor or stator blade 


rows, computations at that refered flow rate were stopped and the Do loop 
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continued. The turbine was found to choke before a refered flow rate of 
4.0 was reached. Of course the point of choking will vary from one turbine 
design to another depending upon the dimensions of the blade rows and the 
losses. A sufficiently high upper limit should be chosen for the Do loop 
so that choking would occur prior to completion of the Do loop. 

The flow function for each blade row was computed and compared with 
the maximum value of the flow function corresponding to the critical pres- 
Sure ratio. If the maximum was exceeded, the name and number of the blade 
row was printed out and calculations at a new refered rpm was initiated 
If the maximum value of the flow function was not exceeded in any blade row, 
the calculation of the performance parameters was completed and the answers 
printed. 

C. Subroutines for Stator and Rotor 

Frequency occuring constants and exponents used through out the 

program were computed in function sub-programs. All other computations were 
made in subroutines Stator, Rotor, Diffusor, and Ratio. As shown by the 
flow charts, the form of the subroutines for the stator and rotor are quite 
similiar. The inlet flow angle was utilized to obtain the loss coefficients 
fora particular blade row. Since the loss coefficients were picked from 
the loss curves at ten degree intervals and presented as a one-dimensional 
array, interpolation for intermediate values of the inlet flow angle was 
necessary. This interpolation was accomplished by subtracting the value of 
the inlet flow angle from 70° and dividing by the 10° interval. The quotient 
must be added to 1.0. The angles must be expressed in radians for all com- 
putor calculations. A change from floating point to fixed point arithmetic 
caused truncation of the result to a whole number equal to or greater than 
1.0. By obtaining the difference between the floating point and fixed point 
values, linear interpolation between the closest given values of the loss 
coefficients could be made. This method of interpolation was believed to be 
sufficiently accurate over the range of angles for which loss coefficients 
were calculated, since the percent error of the loss coefficients was not 
known. 

D. Subroutine for Determination of Pressure Ratio 

In order to obtain the pressure ratio across a blade row a separate 


and rather complicated system of tests and calculations was necessary. A 
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subroutine called Ratio was formed. The value of the polytropic exponent 
was calculated and used to find the critical pressure ratio for a given 
blade row. The critical pressure ratio was substituted in Equation 5 and 
the value for the maximum flow function existing at the critical pressure 
ratio was determined. The calculated value of the flow function was com- 
pared with the maximum value and if the maximum was exceeded, the blade 
row was choked and a return statement would transfer control back to the 
main program. 

In order to find the pressure ratio corresponding to a given value 
of the flow function obtained from Equation 5, an approximate value of | 
the pressure ratio was computed using Equation 6. The approximate value 
was tested to determine whether it was greater or less than the critical 
pressure ratio. If the approximate pressure ratio was greater than the 
critical, .05 was subtracted from the initial value and the new approxi- 
mate value, which would be less than critical, was used. If the approxi- 
mate value was not greater than the critical value, the approximation was 
increased or decreased in steps of .0001 as necessary, and a trial value 
of the flow function calculated and tested at each step. As soon as the 
known value of the flow function was bracketed, the last value of the 
approximate pressure ratio was considered sufficiently accurate to use 
as the pressure ratio corresponding to the value of the flow function. 
In the case of manual calculations, the pressure ratio could be obtained 
from Table III by making a two way interpolation. The use of logarithms 
was required in order to achieve the necessary accuracy. Since only sub- 
sonic flow has been considered, all pressure ratios will be less than the 
critical pressure ratio. 

A flow chart of Subroutine Diffu (Diffusor) is not presented since no 
control transitions or conditional statements were involved. Only a 


Straight forward step by step solution of the equations is required. 


VIII. Turbine Analysis 
A. Preliminary Analysis based upon Design Drawings 
In order to demonstrate the application of the method and the 
usefulness of the computor program, 15 different computor runs were made. 
A list of the runs is given in Table VIII, Appendix V, which shows the re- 


fered rpm, specific heat ratio, gas constant, rotor tip clearance, and 
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blade flow angles of each of the runs. Print outs of the results of the runs 
are given in Appendix \V. 

The turbine had been designed to run at 18,000 rpm and 1200 °F using 
Nitrogen gas as the working fluid. Since the test data presented in Table 
II and III were not available at the time of completion of the programming 
of the equations for computor solution, a series of computor runs were made 
using flow areas and rotor tip clearances based upon design data and the 
blade drawings. 

Exactly how much the rotor tip clearance and corresponding flow areas 
change with changes in temperature depends upon the design and the materials, 
and cannot be determined except by extensive testing. It was assumed that 
the high temperature of 1200 “F would cause both rotor tip clearances to 
be reduced, since some thermal expansion of the rotor blades and walls would 
take place. The turbine design indicated that the reduction of the clearance 
of the first rotor would be greater than that of the second, therefore the 
clearance of each rotor blade row was assumed to be the same at high tempera- 
PuEe Operating conditions. , 

1. Test Runs at Design Refered RPM 

The first four runs were made at design refered rpm using Nitrogen 
as the working fluid. The clearances and areas measured from the drawings 
were used for Run #1. The rotor tip clearances of .033 and .021 were reduced 
to .015, .010, and .005 for runs number 2, 3, and 4 respectively. The re- 
duction of the tip clearances caused a corresponding reduction in the flow 
areas as shown in Table VIII. 

A plot of refered flow rate versus pressure ratio was made, Fig. 33, 
using the results of these runs. The plot showed that the pressure ratio 
required for a given flow rate increases as the clearances decrease. 

2. Development of Maps and Indicated Turbine Performance 

Based upon the very limited information available, a rotor tip 
clearance of .015 was believed to be the best approximation of the actual 
clearance that would occur during operation. This clearance was assumed to 
exist over the range of refered rpm from 240 to 540. For Run #5 the com- 
putor was programed to compute the performance parameters over the complete 
range of refered rpm at intervals of 10. The refered flow rate covered the 


range from 2.0 to choking in steps of .1. Since a relatively large increase 


a] 








in pressure ratio is required to produce a small increase in refered flow 
rate as the pressure ratio corresponding to choked conditions is approached, 
an additional run, #6, was made over the same range of rpm, but over a re- 
duced range of refered flow rate starting at 3.7 and increasing in steps 

of .01. The performance parameters corresponding to flows near choking were 
more accurately determined yet the data output and computor time involved 
was not excessive. 

The complete set of performance maps, Fig. 26 through 32, were drawn 
using the results of runs #5 and #6. Although the data input to the com- 
putor was based solely upon the design information and an assumed rotor tip 
clearance, a general overall prediction of turbine performance can be made. 
This estimation could be refined when actual measurements of the blades and 
clearances were made. 

The three most important graphs of this set are Fig. 26, 27, and 28. 
From these maps all of the parameters which are needed to completely define 
turbine performance can be obtained if any two parameters are known or 
assumed. Turbine efficiency and power coefficient corresponding to given 
values of pressure ratio and refered rpm are presented in Fig. 26. The re- 
fered flow rate can be obtained from Fig. 27. The velocity ratio corres- 
ponding to the efficiency and pressure ratio of the operating conditions can 
be determined from Fig. 28. 

In order to draw the curves of constant pressure ratio and efficiency 
in Fig. 26, plots of power coefficient versus turbine efficiency, Fig. 29, 
and power coefficient versus pressure ratio, Fig. 30, were made. The data 
sheets for run #6 were also consulted in order to determine accurately the 
curves of constant efficiency for .780 and .784. 

From Fig. 26 it can be seen that turbine efficiency is a maximum of 
./84 at the design refered rpm of 441.8 and a pressure ratio of 2.69. The 
power coefficient corresponding to these condtions is .189. The effects of 
changes in either the pressure ratio or the refered rpm upon turbine effici- 
ency can be clearly seen on this map. 

A refered flow rate of 3.72 was obtained from Fig. 27 for the design 
refered rpm of 441.8 and the pressure ratio of 2.69. Fig. 27 shows that 
the pressure ratio changes for a given refered flow rate with changes in 


the refered rpm. This effect upon pressure ratio is more pronounced at the 
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lower values of rpm. The refered flow rate and pressure ratio at which 
choking occurs for a given refered rpm is shown by the point of termination 
of the upper ends of the curves. 

Fig. 28 shows the velocity ratio corresponding to pressure ratio and 
turbine efficiency. This map was developed from plots of velocity ratio 
versus pressure ratio, Fig. 31, and turbine efficiency versus pressure ratio, 
Fig. 32. The curves of constant pressure ratio in the figure show that the 
Operation of the turbine will be restricted to a fairly narrow range with 
the maximum efficiency occuring at a velocity ratio of .445 for a pressure 
ratio of 2.69. This plot also shows that the efficiency varies very little 
until pressure ratios too low to be used in normal operation are reached. 

B. Performance of Turbine with Redesigned Nezzle Blades 

The manufacturer had indicated that the nozzle blades would be 
redesigned so that the flow would enter the nozzle blade row at zero inci- 
dence. Such a redesign would reduce the losses in the nozzle which were 
at present excessive due mainly to the large negative angle of incidence. 
The flow areas of the Nozzle, Rotor I, and Rotor II were to be reduced to 
8.28, 9.38, and 14.20 sq. in. The planned modification would reduce the 
loss coefficients (_ and C from .2050 and .2475 to .0828 and .1364 re- 
spectively. 

Run #7 was made using the areas and loss coefficients for the modi- 
fication so that a comparison between the turbine performance of the ori- 
ginal design and that of the redesign could be made A plot of refered flow 
rate versus pressure ratio for this run was included in Fig. 33. Froma 
comparison of the printout of Runs #7 and #3, it can be seen that turbine 
performance would be improved by the redesign. The maximum efficiency would 
be increased from .791 to .819. 

C. Attempted Correlation of Test Data and Program Predicted Performance 

When the open cycle test data presented in Tables II and III were 
received from the manufacturer, several computor runs were made at the same 
refered rpm as the tests in an attempt to correlate the theoretical results 
and the test data. There was insufficient time for a complete comparison 
to be made. A complete set of performance maps was not drawn in each case, 


only a plot of refered flow rate versus pressure ratio was made 
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The test data presented in Table II was obtained by the manufacturer 
from open cycle tests using air as the working fluid and methyl alcohol as 
a fuel for in-line combustion. The turbine inlet temperature and rpm were 
slightly lower than the design values and resulted in a refered rpm of 
420.9. The pressure ratio across the turbine was 2.347. The mass flow 
rate was taken as the sum of the fuel and air flow rates and was 3.44 lbs/ 
sec. The corresponding refered flow rate was 3.925. 

In order to determine the specific heat ratio and gas constant of the 
flow, the ccmbustion gases were considered to be the products of the com- 
plete combustion of methyl alcohol and air. Sample calculations of the 
Specific heat ratio and the gas constant are presented in Appendix I. 

The manufacturer had measured the actual minimum flow areas of the 
blades assuming that a rotor tip clearance of .020 existed when the tur- 
bine was operating. The measured areas were different from those obtained 
from the original drawings. The measured areag are listed in Table II. 

Runs #8, 9, and 10 were made using loss coefficients and flow areas 
corresponding to rotor tip clearances of .010, .015, and .020 in. re- 
spectively. A plot of refered flow rate versus pressure ratio, Fig. 34, 
was made from which it could be seen that the measured flow rates for a 
given pressure ratio were greater than the theoretical at all three values 
of tip clearance, and that either larger clearances or flow areas existed 
in the operating turbine. 

Runs #11, 12, and 13 were made using a refered rpm of 407.4 which 
corresponded to the temperature and speed of Test II. The refered flow 
rate determined from the test data was 3.995. Due to the simularity of the 
conditions for Test I and II the curves of refered flow rate versus pressure 
ratio shown in Fig. 34 plotted very close to those for the higher rpm of 
209 . 

Additional test data was received from the manufacturer as the open 
cycle tests were completed. The refered flow rate and pressure ratio of 
each test was calculated and plotted in Fig. 34 resulting in a scatter of 
test points through which a single average curve was drawn for comparative 
purposes. The curve indicated that greater flow rates for given values of 
pressure ratio actually occured than was indicated by the computor program 
results. The percentage difference was not as great as the expanded scale 


of Fig. 34 would indicate. 
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Test data for two tests which were conducted at low inlet temperatures 
was also included. The first of these test, Test III, was conducted at an 
inlet temperature of 715.5 °R. The low temperature resulted in a tempera- 
ture drop through the turbine of only 150 oR Small errors in measurement 
of total temperature under such conditions can result in large percentage 
errors. The test was made at 14,000 rpm or a refered rpm of 523.4. Com- 
putor run #14 was made at this refered rpm using areas ratioed down from 
the measured values to values corresponding to a rotor tip clearance of 
.015. A curve of refered flow rate versus pressure ratio was drawn, Fig. 
34. Run #15 was made using the refered rpm of 594.3 at which Test IV was 
conducted. The results of the computor run was plotted, Fig. 34. In both 
cases the flow rate obtained by the actual tests was slightly greater than 
the theoretical value calculated by the computor. The efficiency calculated 
from the actual tests was much greater however, 84% compared to 77% for the 
theoretical computations. 

During the final reading of this paper it was discovered that the sign 
of the incidence angles of the flow into Rotor I and Rotor II was not in 
accordance with the sign convention adopted, Fig. 2. As a result the lass 
coefficients presented in Tables IV and VI, and the corresponding graphs, 
Fig. 18 and 20, are in error. The error involved will not cause a signifi- 
cant change in the loss coefficients or performance of the turbine under 
normal operating conditions. At reduced flow rates, where the incidence 
angle of the flow entering a rotor blade row is considerably larger than the 
blade angle, the loss coefficients will be less than those presented in 
Table V and VII. 

Vavra continued the investigation of the performance of this turbine 
using the same basic method. Loss coefficients for the rotors were cal- 
culated with the sign of the incidence angles of the flow into the rotors 
taken in accordance with the sign convention presented in Fig 2. Since 
the Reynolds number of the flow corresponding to design conditions is approxi- 
mately 7 x oO and is considerably greater than 2 x 10° for which the data 
in Ref. 6 applies, the profile loss coefficients were corrected for Reynolds 


number effects using the empirical relation suggested in Ref. 6 


Cp = yay xp (for R, 3 2x 10°)| 


ot 









The values of ( p were reduced by 22.2%. 


Vavra developed a computor program independently and made a run using 
the same assumed rotor tip clearances, flow areas, and refered speed as 
Run 2. The results of this run are presented in Table VIII. A maximum 
turine efficiency of 80.6% was obtained for a pressure ratio of 2.56 and 
a refered flow rate of 3.80 Considering the reduction in profile loss 
coefficients, the efficiency compares favorably with the value of 78.4%, 
for a pressure ratio of 2.69 and a refered flow rate of 3.72, obtained 


from Run 2. 


XI. Conclusions 

This method of turbine performance analysis will provide an accurate 
and rapid means of determining the performance of a subsonic, axial flow, 
multistage turbine providing the actual measured flow areas and rotor tip 
clearances existing during operation at high temperatures are known. A 
reduction in clearance will cause a corresponding reduction in flow area 
and flow rate for a given pressure ratio, and increase the efficiency. 

All the dimensionless parameters needed to completely define turbine 
performance can be obtained from the turbine performance maps. These maps 
show the effects of changes of one or more parameters upon the others. The 
plot of refered flow rate versus pressure ratio shows that a relatively 
large increase in pressure ratio is required to produce a small increase in 
refered flow rate at pressure ratios close to critical. 

As a result of the analysis of the performance of the two stage turbine 
investigated it can be concluded that the turbine operates at a maximum 
efficiency of 78.4% when running at the design refered rpm and a pressure 
ratio of 2.69. The power coefficient corresponding to these conditions is 
-189. The operation of the turbine will be restricted to a fairly narrow 
range of velocity ratio with maximum efficiency occuring at a velocity ratio 
of .445. The efficiency at this velocity ratio varies very little until 
pressure ratios too low to be used in normal operation are reached. 

The redesign of the nozzle blade row, as proposedby the manufacturer, 
would reduce the pressure losses and increase the efficiency approximately 


3% to a maximum efficiency of 81.9%. 
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Correlation between the test data and program results was as good 
as could be expected considering the limited amount of test data available. 
In all cases the measured flow rates for a given pressure ratio were greater 
than the theoretical values determined by the computor programs, however © 
the percentage differances were not excessive and closer correlation should 


be possible with additional computor runs 
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Blade tow end Gas Angle Geometry 


¢ chord -- straight line connecting the ends of the camber line 
t maximum tnickness of blade 

S$ spacing -- blade pitch 

a blade opening or throat -- height of minimum flow area 

t, trailing edge thickness 

&<{4 inflow angle 


of discharge angle 


qn 
‘\ 


blade angle at inlet 
({ blade angle at trailing edge 


i incidence angle 


A 


designates stator blade angles 


Ly 


designates rotor blade angles 
Sign Convention: (1) Angles are positive where velocity vectors have 


components in the direction of rotor motion. 
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\2) Incidence angles are positive when the def- 


lection angle is greater than for a flow 
envering at the blade angle. 
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Fic. }/4. Positive stalling incidences of cascades of turbine blades. 
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Fic. ! 5, Variation of profile loss with incidence for typical turbine blading. 
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Fic. 1/6, Secondary losses in turbine blade rows. 
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TABLE I 


BLADE DIM=NSIUNS AND ANGLES 


Inlet Diameter 

paeenanes Diameter 

Mean Average Diameter 
Mean Blade Height 

Radial Clearance 

Maximum Profile Thickness 
trailing Edge Thickness 
Number of Blades 


Throat'a" 


Throat Area sa. 


Chord 
Spacing at Mean Diameter 
Inlet Blade Angle 


Outlet Blade Angle 


i 
in. 
aba 
ait: 
in. 
tye 


ie 


vnc 
1vic 
ire 


in. 


Ges. 


deg. 


NOZZLE 


12.0 
tO 
1276 
660 
O 

. LOO 
O21 
79 
a9 
Saou 
. 660 
207 
ee) 
70.0 


nOTOR 1 


Vo 
Ley 
Laee> 


s{O5 


proU 
Ele 
100 
4865 
1022 


Cee 


te 
eu 
oO 
ACHES, 
0 
5098 
020 
io 
<oO> 
jy 
eer 
5 5de7 
ah 1, 


W2eO 


Owen, Li 


pa 
ee 
le 
1.004 
Oe 

nOOT 

~034 

83 

. 1646 
ee 
2765 

-498 
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OPEN CYCL® Test DATA 


Working Fiuid--Air; Fuel for In-line Combustion--Methyl Alcohol 
See TES ee 
Turbine inlet Pressure psla. ahs 250 
Turbine Uutlet Pressure psia. 4.7 Lae 
Turbine Inlet Temperature On 1090 17ey7 
Turbine Outlet Temperature = 825 917 
Turbine Air Flow lb/sec aon 3.30 
Turbine Fuel Flow 1b/sec 0901 1027 
Compressor inlet Pressure psia. fie, rae 
Compresco. OULLeE Pressure psia. 36.5 a7 oO 
Compressor inlet Temperature OF 60 64 
Com -.ssor Outlet Temperature OF 281 284 
Compressor Air Flow lb/sec 3.559 3.490 
Turbine Speed rpm too 10 16,620 


blow Areas as measured by the manufacturer assuming a radial tip 


clearance of .C20 inches; 


NOZZLE Ons. SQ. Ing 
ROTCR I 10.55 sd.cime 
STATOR Hi./2 sQGluoe 
ROTOR We 15 LO Soames 





TABLE IIT 


OPEN CYCLE TES? DATA 


Working Fluid--Air 


SDS Ma BS TiS ay 
Turbine Overall Pressure Ratio 2eOay 2 soe 
Turbine Inlet Temperature oe Tia 1i5ee 
Refered Flow Rate on a. age 
Refered RPM D234 Soa 3 
Turbine Efficiency % Sloe 84, .00 
Additional Test Points 

Pressure ratio Ref. Flow Rate Pressure Ratio Ref. ‘ow Rate 

1.28 Dwee/ ae 3.22 

trons 3.08 1.28 fe) 

oo 3.46 1.96 3.63 

ae alo | Zon See 

aa pee 2a 0 3.875 

1.47 S250 Bea? 4.00 

LeLe, 3.38 

53 Belk | 

oo 2205 | 

1.65 Bec | 

Oo Eni 2 

oO E207 

1.80 3.96 
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APPENDIX I 


Development of Equations 











I. Develooment of expression for the Polytropic Exponent 
in terms of Ye 
The efficzemey of the flow protess from the amtetsto tne minimum 
area of a blade row is equivalent to {1 - ae and can be expressed 
in terms of the temperature ratio éT/dT; . as 
YY = at/at,, = 1 - 5 


In general 


Aoi 

(T- aT;,)/T = ((p - dp)/p) 
a) 
1-d7,,/1 = Ga ep/p) 


Expressing the right side ot tne equation in series form (1 -€) 


L- aT, = 1-(f%-1)/p ap/p 


or aT; /T = (ia) aeiece © 
Since 
arf= dT.) 
then ar/t = Y) CS =) 
Intergrating An} 
InT = Inp’#F ¢ Constant 
sent 
thus T = pf oe 


In order to express this relation in the form of py" = Constant. fom 
a polytropic process 
QAY- Vie = GsaGa- W/m = (n - 1)/n 
or i/n = (1 Sie) 4 )/r 


and n 


PF (8-1) #2) 
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II. Derivation of the ecuation for the Approximate Pressure Ratio 
across @ Blatc wo: 


The Flow Function was defined as 





mae rettoup/P., is equivalent to (Pi =Ap)/Po or i - Ap/P.. 


Expressing the exponential terms in a binomial series expansicn 


(p/p,)*/” = i+ a AR uy (2/n -1) (“3 y* Le ee eas 
ho. n a 
(p/P we f i)/n 
O 


_ ri ae n 71 /Ap Z 
= to aes) 


Substituting the series expressions in the Flow Function Equation 


and reducing 





mo 7 a } = Ap \* 
Flot Pte + 302909) 


oO n 


Solving tor Ap/? , and reducing 


. >> 
\| (ea) <— 
— f= 7) 


Ap/P, - $B 
3 





so for a lst approximation P./p Cen Se expressed as 





Eee = L = n | oe oe o* i) 
: I - apr, y}i- Bip. oP Gi i f 


O9 





III. Sample calculation of the specific heat ratio and the gas constant 
of the combustion g“ses oroduced by the combustion of Methyl 


Aleehed: and Air, 


For Test I: Fuel Flow = .0901 lbs/sec 


Air Flow = 3.35 lbs/sec 
MW.OW. of CH) OH =o. 
Mie We Cl Air = 25.97 


A/F = 3.35/.0901 = 37.2 lbs Air/lb Fuel = 42.35 mols Air/mol Fuel 
Basis: 1 mol Fuel 


CH,OH # 8.9 Oo # 33.45 No ---- COy f 2.5 HoO f 8.15 0, # 33.45 No 


20.8 Oxygen atoms available 4.5 Oxygen atoms required 
Therefore 20.8 - 4.5 = 16.3 Oxygen atoms excess 

or 362% excess Air 

By entering Table 5 of Ref. 5 with the average temperature of 1420 OR 
a specific heat ratio of 1.344 can be obtained. 

The average mol. wt. of the products of combustion was found to be 
28.7. By dividing tnis average into the universal gas constant (mR = 


1545) a value of R of 53.9 was obtained, 
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APPENDIX I1 


Loss Coefficients 


Sample Calculations & Tables 


a 








BOUATIONS AND CALCULATIONS JNVOLVsD IR LOSS COMPUTATIONS FOR NOZZL:s 
From Blade Profiles: s = .509; c = .600; t = .100; t,= .021; 

a= .1y0; X4"2 70.09; = -h2.5; 

sje .7713 t/e.=/2oe t/a Se eo 


mom Tunbane Drawing: ee .ood ain, 


—__ =P = =P aw ee geen mae ee Ragen or nn a ere Po gece Me gem am 


= 61,89 (<4") 


ne 
BG 0) = Yo (a gt QO) Zc; |: ‘py yy “pegs “of (wg) ( 


=f 030 ee 5)" | .126 7 -030/¢ ( Igie)" = ae = = cer 


easyer ay a 72> a 61,.8/.995 = 65 .0° 
Xo ek i(s/e = .75) = ~42-5/09 = ~.0545  gis/ex 75) = 20.0 


Ai.= -1.2; a= 20.0 f (-1.2) = 18.8° 


\f 


Ss 
iy de = -42.5/18.8 = -2.26; tp/Lo(is Q) ~ Bel. toi ey ame) = Leo 
Bere fp/) 4) = 1995/1 eee 

fee 162/152 ) = 1.2/(l eee cio 

£5" /az $/\1.2 #3) = .1670/(1.2 # .1070) = .122 


ie St eo ~19* - Yxil- 897°) = 2050 
Aa = TD nycosXy” = 3.14% Poe ee core (LI. DS 
A, = TT Doh pcos Ky = 3. i wee oom = 1130 


j 2 
: A) _ » Liye CT. L67 
A = r| 7 ay - x Le 3-878) | = 8 [ J 


Saee07O° trom Fig. 
Kn = tan-+ (tan &%, ¢ tan & ee ~ tant (O # 2.13)/2 = L6.8° 


2 
L A(cos 1/cose_) (tan&, - tan, ) 


(t 


eSseel2 
hx .0076 x (Bey geer © — 2.13) = / ¢03en 


if 


‘. pete ey Y_) =. 0304 / ge Ooo ee 0351. 
ies = Bi te/2)(3,/2€ /a) See eee 67 /.122 = 0154 
Srotar® 9, 435, 4 V2 = +1070 4 20551 4 044 = 2475 


2 





BOUNELONS ASD CRICULETI Geer VOLVO) in LOSS COWRUTATIONS FOR ROTOR I 
Prom Blade Profiles: s = «4665; ¢ = .760; Usamye 7c, 2.) — .030; 

dea 605 As" = -70.83 Ay a2 

sfc = .639; t/c = .17833 to/a = 2275 
From Turbine Drawing: K(cold) = -0335 h = 5705 


—_ ews wm sen wu i ss owum ise ose ism ose ose ise ia oss «es_# am == ss ss = CSP SEP Ss SP SP Ce ee OleOle 


[So = cos7 a/(s - te/cosfs ) = cos~l .160/(.4865 - .036/cos-70.8) = 


if 


-64,..9° 8% 
(2; 


| 
mpta = 0) = ftpisi"s 0) #(3 Veteay? tp ato] (ele) 
= }.034 + ( 1008 2° [227 - 031] p (1202) Fates = 050 X 690 ate 


6L,.9 
Bo(s/o = 75) = *2/1.035 = ~64.9/1.035 = -62.7° 
/3\"fB2(s/c = 75) = 10.2/-62.7 = -.1033 is(s/e = .75) = 35-0° 
Ai 3255 ig = 35.0° # 3.5° = 38.50 
Bins ie) ee Yn) 80 ee eae, 
R= 2/2! 5 5) = 12/ Cee cca 
Ser 09 (1-97) = 9 x (L = 29742) = 046 
£8" fa = Sp/\1.2 £¥,) = .0324/(1.2 # .0324) = .0202 


Ay = (1Dyh\cos 4) 20321), x 1258 oe ete oseloee Bao.) 


Ao = {Donocos 5 = Bly x 129 ee (ere ser oy =e 
nee x [soe ie y= F fe .9/26.1)* ear] = fee | 
Hl # (ID/OD (12.8/12.9). 


mWoS 1rem Nig 7. 
Equations for Y, ana’Y, nold between -1.5 and 1.0 i/i.: (See Ref.-6) 
Yeas hA(cos*4 o/cos f3,.)(tan/ ; ~ tan/7,) - - 

= x 0008 x (.180/.690) x ( 0 - (-2.41))"= 10326 


2 k/h (cos ,/cos 4m) (ten4A, - tenfio)* = 


Yi 
2x .033/.705 (.180/.690) ( O- (-2.11))* = c1ey 


where ae - ten-L ((tan 74, f tan #5)/2) = tan7+(-1.055) = -4o.4 


ee 











Yeu /tb #1) = .O3L0/(1 # .0316) 
Sk ca Y/(L ¢ Beye 1090/{1 #4 .1090) = .0990 


Fee *Blte/alp/2S/a) = «3 x .225 x .0324/.0262= .083 
Nee ie 7, 4 ‘. £ ,,. fi eo = »O37k fs ,0306 f 20990 f 00835 = - 2460 


7h 
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APPENDIX III 
Basic Fortran Program 


Table of Fortran Names and Symbols 


Flow Charts 


78 




















ALPHRO OC 3) 


a *e(2,4) 
AS Ae(1,3) 
BETAL P13) 
BETO Poo 4) 
cl Cy 

C2 C. 

03 C.3 

CP Cy 
DIFF 

DLFFU 


DUT Din 1,3) 
DMO Da(2,4) 
vTD el / lo 
Mitso, AT. ./T, 
LTO Ar], 


DTH ATp/T, 


UTS ATS/t, 


FORTRAN WAMES. KE ulv’ Lin? SYVBOLS. AND DEFINITIONS 






ntermediate step in calculation of the Approx- 
imate pressure ratio acress a blade row. 


Angle of absolute velocity et inlet to Nersie or 


Diarore 


Angle of absolute velocity av exit from Nozzle 
or Stator. 
Minimum tlow area between rotor blades. . 
Minimum flow area between stator blades. 
Angle ot relative velocity at inlet of Rotor, 
Angle of relative velocity at exit from tovor, 
Uonstant; See Sample Calculations 
Constent; See Sample Calculations 
Constant; See Sampie Calculations 
Speer: gc eee at Constant pressure. 
Ditierenice 
Diffusor 
Mean diameter at inlet to rotor blade row. 
Mean dlameter at exit from rotor blade row. 
Issntropic temperature drop through the cifrusor. 
isentropic temperature drop across a blade row, 
Actual temperature drop across a blade row, 

; 
Temperarure ris ue LO the kenetic energies of 
the relative and perpheral velocities pertaining 
to tne rotor. 


Yemperature rise due to tne kenetic energy of 


the apsolute velocity of the flow leaving the rotor. 


, 





DTT ATs s7/lo Isentropic temperature drop through the turbine. 
a> s M2 ee : e + 
LUG 1/\P/p) intermealate step in the calculation of an 
Approximate Flow runction. 
DUGH i/r/p)” lnto¥medlats etep if pc Galeuldtlon ef ths 
critical value of the Flow Function, 
my ; . f A 2 
DUN 1/(P/p) intermediate step in the calculation of an 
Approximate Flow runction. 
my) . 
UUNM A/D). Intermediate step in the calculation of the 


critical value ot tne Flow Function, 


DTW At, /% Actual temperature drop across the turbine blade 
rows. 

aN n rolytropic exponent 

wTAD No Diffusor eificiency for ideal axial exit er tne 
LLOW, 

ETADA Ns Actual aitfusor efficiency. 

aay .. Overall turpine efficiency. 

BXPL 2/n Pressure ratio exponent for polytropic process. 

EXP2 in # 1j)/n_ Pressure ratio exponent for polytropic process. 

EXP 3 m5 PYessueecm@aeto emeonenL LOr ISssntrople processr 

wxPL Mo Pressure ratio exponent tor isentropic process. 

EXP 5 n/(n #1) Pressure ratio exponent tor polytropic process. 

GAM 7 Specific heat ratio (C/Cy). 

TER Flag allowing decision for control branching. 

ICR Flag allowing decision for control branching. 

IFLAG Flag allowing decision for control branching. 

i Pass or stage number. 

Or @ Value of tleow function correspcnding to the 


approximate pressure rstio. 
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OTA 


OTM 


PC 


PR 


PRA 


PRE 


PRHO 


PROE 


PROR 


PROS 


PRRO 


PRSO 


inf 


RRATE 


RRPM 


{H 


(P/p) a 


(P/p)e 


Flow Function for a blade row. 

Maximum value oft the flow function for a blade 
row; choking soccurs, 

Ratie ef the static pressure after a row of blades 
to the turbine total inlet »ressure. 

Horsepower in coefficient forn. 

Ratio ot total pressure at inlet to static press- 
ure at exit for a blade row. 

&eproximate total to static pressure ratio across 
a blade row, 

Critical pressure across a olade row. 

Ratio ot total pressure at exit from diffusor to 
static pressure ahead of the dirfusor. 

Reciprocal of the overall turbine pressure ratio. 
Overall turbine pressure ratio. 

Ratio of total pressure after a rotor blaae row 

to the total pressure at iniet to the turbine. 
Ratio of total pressure after a stator blade row 
to the total pressure at inlet to the turbine. 
Ratio of the total pressure to the static pressure 
at exit from a rotor blede row. 

Ratio of the total paeecure to the stetic pressure 
av exit trom a ststor blade row. 

Gas constant;(lo45/ifolecular weight) 

Refered Flow Rate 


Refered Speed 


Noncimensional total temperature atter a blade row, 


SL 





"TS 


ate 


TS 


Ul 


UO 


VRATIO 


yU 


WU 
TS) 


ZEUS 


ee 
Tr/T 
T/T. 
Bir 3 IT, 
Yo, 


Uf le 


Vp/}T, 


oe 


Unya./Co 


Ratigot the temperature at inlet to the temperature 
at exit of a blade row for lsentropic conditions. 
Eouivalent temperature at inlet to a rotor blade 
row . 

Hquivalent temperature at inlet to a stator blade 
row. 

Perpheral speed at tne mean rsdaius and rotor 
inlet. 

Perpheral speed at the mean radius and rotor 
exit. 

Absolute velocity of rlow. 

Absolute velocity of tiow at discharge trom the 
dit iusom, 

Meridional component of the absolute velocity. 
Velocity ratio tor the turbine. 

Theoretical velocity for isentropic expansion frrom 
Stagnetion pressure at turbine inlet to static 
Pressure ev (GMetpsoree eciarsc. 

Perpheral component oz the absolute velocity. 
telative velocity of the flow. 

Perpheral commer saemeneure relative velocity, 
Total loss coefficient for a-row of blades. 
Expansion loss trom inlet to the throat of a 


blade row. 
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Heaa Inout 






Write Heading J 
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| Do RRATS .1 ~- Choke (.1) | 


| 
Initialize } | | End 
TS; PROR: 7E: 7: ALPHAL / No 
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Flow Choked 
in Rotor No. 
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SUBLOUT UI] Re 


Stace 






| Dele = Zeno 


| Calculate EN | 


Calculate PRC | 


Y¥ 


Calculate Of 





¥ 


Yes se ade ory 
| 


NO 
| 


Calculate PRA | 





tes PRA » PRC 







Subtract .05 | 
Peon - RA 







Increase PRA 


. 


Decrease PRA—>< Yes 
by ,0002 











f Equal 
Xecalculate A Recalcusave 
nt 
Set PR = PRA \ 
Yes | ~ No. OT > OTA 
No Yes 


| : { 
<p Return Set PR = PRA 


89. 





ECUATICNS #OR CRICUL: ‘TIONS MADE IN SUBROUTINE RATIO 


CUS erate v 
SUF- TIT 


(2) (P/P)eritical = in Z 1.)/2{ n(n - 1.) 


Oye? fan yaa. f(y - (33) , mAl. | 


—_- LS P/p crit. P/e/erit. 
(4) (P/p); = ie 
Approx. a yn ae ye) Wi See Zz) 
3 } wpin=- 1.) 45 
) i, 
(5) =e = pd (2.,-) /n - (dn \ # 1.)/n 
“ -~ 1, Vee Approx. 7 hopron 
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APPENDIX IV 


Sample Calculations 
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SAMPLE PERSORYANC: CALCULATLONS USING VALUES 
YP Yok OPEAATING rARAMMTERS WeIChee ee 


MAXIMUM EFFICL“NCY FOR THa DESIGN nSabveRe) SP#ED 


Rigg. 10 ft/ OR ; To = 1200 °F = 1660 OR 
P= 1.36 ; No= 18,000 rpm 
wT =a) idm YOR Rotor tip Clearances = .005 
ie sec psia 
N= 4bl.g rp2 
z7 oR 
Ceol  1.36-1. 2 eee 
= 1.36 
mo. _  _ _1.36  _ 3.77778 





pel 1.36-1. 





Gime jE. 755.16 = 1.30936 
ge | S2en 7k 


Co . fom . V2 x 1.36 x 32.174 x 55.16 _ 115.80 
i 1,362 


ee ex l0+ 1. x 104 a ote 


— 


2¢erf/(y- Ll). ~ 2 x 32.17h x 55.16 x 1.36(1.26 - 1) 


TT = .0043633 
720 


Cross~sectional Area of Diffusor ~ 92.3 sca. in. 


Average mean Diameter of Blades - 13.0 in. 


SA OMENS NE NAN NA tg Ne ne NG SANE BY ee NE Ne NZ te SS ee RS Ne OE NS NG SS NS gy tke Re NP Ns Ne WSs as Rs we 4 Ne NP bs $e ts ws Ne ae ye Nes 
SE IEE TE BOOS Te BODE TOTE DY OC EOTE Fo TITS DOPE EGOS ON IE GE POTS PEO TE POTTS ORCI COSI CTIA C IETS OTC TOTO IN ICD TEBE 
re 
NO77Z UL 
4 sd 
Ll Ea 


x 
Ye = -0765 . Jn = .2475 loss coefficients for Nozzle. 
J 


Fy — wi To YR/g Hobe — 3.9 x 1.30936 x¥L _ .51529 
Po Ae} Po Po 9.91 - - 


bows .22002 

Pl “later Serrectec to .2050 
Pi, = «81966 

Pee eave 2002 


Jee 





To .. ( Bo 0) = (1.22002) *20471 _ 1.05405 
Pl 





Ate Ge To /T sade 7. 2 Orb0o = de, ORE 
ATy _ ATigy (lL. - Fy) = -05128 (1, - .2475) = .03859 
EtG 

ns a 

Ss Sen =. 18 (603859 = 2.74684 
iY a 

Uy = .0043673 N D, = .CO436033 x 441.8 x 12.8 = 24.67h04 
|Top Yo 
Us = -0043633 N Do = .0043633 x 441.8 x 12.9 — 24.86741 
VTo eyes 
a sing) = 2?.7L68) x sim Gime = 20.76725 

et 
Mm. | a) cose] = 22.7/681°% Cos Gh.6 = 9-605000 
}To ~7Tto 
It5 JT Lo 


Wy Yar" £ Way — \J(9.08500)* 4 (-4.09261)* = 10.51,21 
yo ~] T To 


2} 
ATiy _ C3 fan? # ut i mn x Lom4 


= 71577 [(10.5U21)2 # (2h.867h1)? ~ (24.67h04)") = .00296 


fa = 2 f AT = +9611 £ .00896 = .97037 

oO O Lo 

>. E Z ain /t 0] 7 ps 4 ose | Bee . 0356: 

Pay = (Ea) (22) = 1.03565 x .61966 = .84882 

Oo a Oo 

Ahi = tant fw oat Be | = tan™> |, 09261 = ~.39981 = ele 
ie Vai/ To 9.68500 
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ROTOR 





Oe eie aT = .L7729 tren ess Cureton cr: 1 
oe wit VR/z YTp/To - 3-9 x 1.30936 x V.97037 . .53003 
oe 
Bo = Pri/Po - 284888 . .68275 
Pam Pa1/Po ~=—«- 1. 28332 
264,71 
Ba = (Bal ee hc) = 1.05934 
T2is \Po 
ie 
ATis2 = THI pee Dois = Tp, jp /Toy ee 
QT QT QT To yee eS 
~ .97037 x (1.0593h -1.) = .05436 
1.05704 
At, _ Aligg (1. - Uy) =  -05430(1. - 17729) = .04472 
T9 To 
To = Tp — 472 = -97037 _ .OhK72 = 92565 
sy Te . 
Wo Cy VAT - 115.8 | 04472 = 2h.1,8821, 
eee T9 
a W> sinfy = 2h.48824 sin(-o4.9°) = ~-22.17580 
(to (oe 
ng My cos@, = 2h.4882h cos(-64.9°) =  10.38787 
EO PT, 
> Wao 1 Uo = —22-27580 FF eeesoyt — 2.69161 
bee 1 }Ty 
2 2 (a. = 
No =¥ Yun 4 Yu = | 610-38787) 4 (2.69161)" = 10.73092 
eT. T 
oO 
2 —l: 2 mes 
Oe = €3 Vo 10 9 = .7h57 eee 0092) x 10% = .00859 
i Te 
Ts2 = 22 y 42e2 = 92565 # .00859 = .93424 
i Te i 
i . P : tiahto me [te # 2008591 3-77778 1 03550 
Po T/T 92565 | 
ae (e )(22) — 68775 x 1.03550 — 70099 
Po Po/ Po 
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x, = tant fy ft to | = tan? | 2e6ghel) 2 eee 


2 
a ie/ 10.38787 
STATOR 

we - 06698 S= +13980 from Loess curves for Steter. 

$5 . wit, \R/e \Eso/To = 3.9 x 1.30936 x 1.93h2h _ .60497 
Po  Aex PS/Po Li54 ~ 70099 

Po = 1.43770 
6 
85 = polo = : aoe eee 


T5 - (2a) 2 03770) ae 1.10087 


ATse. _ B59 /ta,6) -1.| — -92h2h | 1.10087 - 4. | _ .08500 
=—=—153 = =32 Asaf Tso/t 4) = ec ae | 


oO To T3535 

AT, = Tiss (1. - C5) = 08560 (1. - .13980) = .07363 

oe ° | 

T; - Isp _ 493 = -9342h - .07363 = 86060 

ae us Te 

Neue Co AB... — 115,80 -07303 ee ots 

—3 = =3 

T to 
Uz = 0043033 N Da Oi C08 3029) Ko oo Loetree 25-29275 
ee iT 

U, = -OOL3633 N Dy) = (.0CEses ieee 5.2 25.4572 
yv _ 4 
ae 7 

Yu3 2 V3 sine, = 31.42218 sin 09° = 29.33517 
fie {2 
Wiz = Wg cc8X, = 31.42218 cos 69° = 11.20062 
Ee VT 

W V3 Uy = 29.33517 - 25.25295 = 4.08222 


fee fT «YTS 

\] 2 2 2 2 
a = V3 - 133 = {/(1L.26062) < Wim. 05222 ) = Lohr 
ae To To 


i: ey, oe ome 
6Tn3 = a1 ya Us" | 40 
cS Me T5 a 


22 





0 ae a zu 
Alay = + ThSTT |(11.97773)? # (25.44572)* = (25425295)? |= 01143 


oO 
Neg = Ty y ATy, = -80060 f 01143 = .87203 
sips i ie 
Pa. i" LAT: ss ge e 4 a 27/7778 05410 
P3 T7Po. £6000 
Pay - (Ex Bs) ~ 1.05110 x 229175) = Booge 
O P37 \Po 





tant 08228 = 24779 = Ioeo- 


Way3/ z = 
| Vaa/ | [1126062 
ROTOR 
ce = .04648 Spx LOL from loss curves for Rotor LI 
mt 


De wit, YR/2 Vr So — 3.9 x 1.30936 x 1.27203 = .60775 
43 51688 





PL, 

ee BglPo = Sie ee 

eo PR3 PL 1.39286 

tpg = (Eas) = (1.39286) 09167 

Te perk 

AT, — Te. (ERo/Tiis) 1.1 = soneos |ogte7 - 2 = 07323 
Tot > TS3 | Stea/ Tae 4 [1.09267 

ey, = A Tyee (le Sp) ee ec OTF) 06137 
Lo To 

Tm . Tao 4, = 387203206137 92 81066 

O To ae 

ie © WAT, = 115.8 } .06137 = 28.68652 

ie i 

a, 2% sin), = 28.68652 sin(-65°) = -25.99877 

ite (To 

act, cosf7, = 28.68652 cos(-65°) = 12.1235» 

iE 5 pt 

Mok = Mo, 4 Ys, = 2599877 FP 25.hk572 = ~.55305 

}To pT }Po 

V iy = V 2 Y 1501926 2 / ae? ; 
ee Ce eee = | (12.12355) # (~.55305)° = 17.13616 
iTS Th) ike 
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Ty = 3 Wy lo _7.14577 x (1e.izele)” x lo _ 0098 
a i 
Ts, = & 4 42, = +1060 ¢ 01098 = .82164 
sus Ty i. 
~~ re) 8 ‘ 
PB, [3s A Sgu/te] _ fae g sguogdl? Te | ayBeare 
—}- = =f == ct S| = - = 
ay 2 LT | “81660 
Ph /P, £27109 x 1.05216 — .39045 
5 = By. a “ ca 
on (Bo) BH 
xX, = tan fy, AT) — tan) [-.55305 [ = 704559 = 2.6° 
| Vind rm |12.12355| 
DIFFUSOR 
My = .70 for Axial Exit 
Yly = Np cos = .70 x cos 2.0° = .69855 
Am, 2). + Top. ce) 2 ene esc 
i pe 
Coe Re = 55D eos pe 20776 
io. 778.11 1.362 hy 
iq WIT, Cy Aly) 1.055eeoes. 9 eameanjemeam lego se lOc rtgoos 
Po 7 
Wp cL wil, R T/T, Wee ee on eect 
IT, Pa Al Bh ro Pee ee 
c 2 a 
A™isy = Cz a Yu, — Wo"! x 1074 
To [To To J 


= .74577 x .69855 x [(12.12355)? - (5.09145)?] 107 = .00631 
a ig ip Th > zg gp or e707 





Pe T/T “81066 
=e e ea (Bs) = ~ 0297 eae te ee, 38212 
Oo L oO 
Poe oe == Lt C= 2 61699 
Py P./Po ~ 38212 
Mp = Le — (Pe eee) 87) 20182 
is 
= Ao, py vie ee 
oe Aine oe 25,62 © eras 


a 





_ 20043673 N Dayg, _l. 


Co Toe Ts ST! To 


= 20043623 SOARS 30 eee 
eS es oe oe 








APPENDIX V 


Programs and Computor Results 











ay Cie 








es 18 0 eae ee 





Ee EE + end Pe wn ieee 
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T° = HovVauy 28° "VT Woot GO T6°6 STO" CT 255 9E°T OUS? =707¢ 
Gt? YT 79° TT 0eO7 TA°G6 S0O* 91° SS f° T B° Ty 
99° TT ar Co °OT TA°6 OTv" aie oS 96°T ells 
L38° 77 79 °TT CueON TA°6 GTO* QT°SS Qf °T Sg 
SsuTMeIg woos | 
“LeeT) 3} sealy ST°ST oe all 81° TT T6°6 TZO**€€0° 9T°SS 9€°T oy T 
Sy aewioy ody Sy _ ay Ny 30URLEST) Y eure WU ° TOU UN 





ae = 














LLIA HIVE 
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SAMPLES OF INPUT DATA 


RUN 1 


Src). < 
eG CSL Se ca 
Mm OACM Se A 
me i) a NN 
» ® % e e 4 
Sa DO OO 
we = NI) 
Peer wo 
Sa NN 
e * e e e @ 
ee, CS 
NIOCOMO 
ener ti NV, 
See ON 
ee @ © »9 @ @ 
oot <> CS 
PaO at th © 
Ott aMo 
ee ae 
e e ® e e e 

Cl Ge 
eee key sae 
ON oe cay) a 
nS ON 
e.U€°8 ee @ @ 
eee mee Leo Ge sg 
Ree A ee) 
a 
Sa 
ee @ @e @ @ 
So NG et oe Ge 
MNO mer MeN 
LT Bh SOUT Ms oO 
e 2 © e e e 
Nees) Go 
a_i 
DSM LA OS 
Se CONOR 
e @ ee @« @ @ 
oe SOm 
ee ee 
PO DOAN 
— @ e @ e e 
moe oo Ge) 
— ISNMAID 
Tormey 
—— COONIM™ 
= e © e e e 
(Nico Go Cone 
OOO 
a. 
eee Ne) a 
e e é e@ @ 
Cla Oe 
Bs GO we IN 
Ris OC car eis) 
a ee NA OS, 
e e @ e e 


_ rNrm 
~ ee nee, ed) esc 


ow MM DOo 
Ja nl 
ee eS oe 
eeee @ 
oh OOOO 
NINO CO ot 
= Oe 
eeeee 
= eee (C= 
WAR = mC O © 
eT NINOUrS 
oe ae ee er 
“— «€ © © @ © @ 
ee Cis | CD CS ae 
ONM OMS. 
CONGR KS . 
OO Oh) Ni 


2 Oe 3 | lac 2 eee 


RUN 5 or 6 


OCoo$wa 
mONNCOW 
MOM CON) a 
O~—OGNIN 
* @ 8 @ @ g 
COogoenO 
Ne er Qi 
CNN OO 
ee SOc 
@ 0@ @ © 8 @ 
OWN 
NV OOOMm O 
ON TOOK MN 
C-OC.--— 
@® © @ @ @ @ 
CODOOCOMD 
MNO FCO 
OP a Ces 
e@?ee0e e080 
ON Ga 
aN) Oo Sr at oe 
IN et 
© e@ 0 @ @ 
IN TL O 
he NYS OAC 
MOW TON 
mK OCON 
ee ee ee ee ee 
ONT ar cn) 
ON SD re at OS 
WAS LU ee he 
2 ee © ¢ @ 
Or Sinica 
DNS~= OOO 
OOM IAN O~ 
@aeee¢e 
Mm MOON” 
Cam ~eyanN 
eOCTomo 
=a NC ColG! 
mm © © @ «6 «CO 
NOW OM NY 
a ON OS 
esOnrCa! 
Oe SOO 
m= e@ e¢ @ @ @ 
NOCWWYO 
g®NOMvoa) 
YO-—-c=—- 
m= NOMA, 
eo @¢ @ @ 
COLO Oni 
O30 SG 
ONIN © at AN 
0 = ee at 
eee 8 @ 
om Co = ee eae, 
MOON 
MSW A 
ee C= OSL) 
é @ *® @ 98 
Cin Tooe 
oh OOOO 
ANIA NIM COU 
© 0 @ @ @ 
—ANOONANDO=— 
NIP NIOWOO 
vel eh silellaitetes) 
SS On 
— © © @ ¢ # eC 
= Ot oom) 
CONM ON Or 
@ON mm st 6 
OS Oe NN 


e204 28-1.13272-1. 13446 


RUN 7 


Cot) Cope 
ee O'S NS Nee 
he YB TIO ND 
Zao 
eee 4 # § 
eo @coacre 
oN OS? aes 
XD VIS OO O 
aac EOE BEN a 
eee @ @ 
Sa eet oy 
NYOCOR- © 
oe Sie 
See — 
ee @ @ @ @e @ 
7 GaN 
KR NDO Sh 
OS aa 
Se 
e e e e e 9 
CNG See 
nee Nigar Ok 
Nae at eel 
ae 
e @ e e e e 
tC SF Oe 
Pe ah a er 
VSO St oe 
© @ © @ @ @ 
GO IN ree 
NOm = Oud 
LAM Su OO 
2 @ @ @® @ ®@ 
Com Goer = = 
a oe oO 
ODM POM 
eee CoO 
® @e @ @ @ 
Reed ices a) 
Na Or 
oo ee 
SVTCcocaq!=— 
— 2 © @ @ @ 
Ch oe 
Raa ee oe 
oO i =f CS 
eee ON Oy 
e @ @ 6 e 
Ne eh ee 
eiNOr 2 a 
ND a CN) 
SN oN 
e e e e e 
Ci 
oe CS 
OM SSL 
renner Ny ONG 
eo «© @ @ @ 
oe OOM 
Nope ONY 
Ne 
e @ @ e @ 
Beso > 
om OOSD 
OI UNG Mh COU 
eeueen oe CN | oF 
© @ © @ 
MNAQOncc= 
= MA~OCOO 
‘SALI SIO CC 
END 
“~~ @ @© © © © oD 
Sno OCN) 
OINM COMM O 
POONA ON c 
ZOnmOrnN— 


RUN 1O or 13 


eee tl 
oe NIC SM coal 
Rae OX eG 
ie CON 
xP 6 RF # # 
Soh 
oy ae Ne) 
gl CW @ es. Wa 
Cee eee CNN) 
eee © @ 
ee aA = ia ie ee, 
MC SOCIUM 
Oe Ni ay ee eanetes 
Ceo Cc. = 
& e e e e e 
Ue Gere 
er oe aN Neem 
eee @ @ 
CONUS ort 
ae 
OLS Sar 
= Co 
e e 2 ® e e 
UIQ Out =o 
eNO Se 
NM LOY Same ie 
eC 
e e 6@ 9 e e 
C00 a ees 
US SSNS La ae 
——SOoN— 
@ 4 e @ e e@ 
CONC On sas 
~~ oo oa 
NOM OT FO 
ae SiGe tx 
e @ e e e we 
Ses Kee ee 
= = 
oT — 
CN Cy COIN. GS 
—— @© @© 8s @ ¢ 
HOvoconm 
et yrON 
Se aes 
— e @ e e e 
Ce eat tee 
CNG Ore 
NS ON SO 
an ON INN) 
@ e e e o 
DA>NOCOM 
eOTOOM 
Qi ia a 
® © @© @ @ 
=O Coy 
NO ees 
I Li eos 
a Ce ON 
eee @ @ 
TNO OOO 
oOWOO 
NAMM ec 
oe @ @ e 3 
NANODATDO— 
Hh MIO Oe 
STINININOOS 
ow 8 6 @ 6 6 @2 
NOMODMNON) 
MNM CONDO O 
CON —DLA 0 
COOFOrANYy= 
~~ @®@ @¢ &® @ @ @ 
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PRCGRAM FOR RUN L 


~ - NJ 
os oe > 
A, Be ers — 
SS if 
—- ~ = 
ee >. o 
oe 2 ON ~ 
ti CoM e — A >< 
Cae oh — NO Gs 
OO wars al ~ — 
~ Ne \! = el ag 
—_— ai) << — a) cd 
loaeed Zl aye ry .— os ee 
CNN Se Sal oe ~ Ge ~ 
—S Ny a ae 
ee ett Ons n~ ~ wo sat ~ 
a mas tl} > Ss ~ et 
as alt IC If al >< < 
ee cy mo UW O x 
fe Or ame Ny SS ~ Oo pee =f 
Agus bat _~ Aw ww es —_ Oni oo i 
me iLi<t ~~ ae eet a w ot i) MN! e 
ON =U) ~ Ile 20 Oe . ho aa ~ 
oad elt, i ae 4 a ran ax A= >) ? — 
Ge ‘i = NN - = e. © ° x 
Oo. f— cCc-_7 7 & ~< ml Pe _ 
ti ov!) ) hae weet eed St i penne Kn a 
mre ~ HII ~ Sen ~ 
a = iN) ame fe NY — Vy ° 
we ett aD nm <{ ‘4 ae <= 2 = 
Zar OL Se re “> a < | 
Ow eS) Wj Sou ae ae, a 
et —— is oS pon © ae ~— < 
meN a7 It ~C oy, <I « C C5 
<< “ClO ed eee ie - Geet < eee ~ ~ 
2 MAL am {| ae Cf e b— cS ~ 
a — 2 ae ™ an ARO = ~ LY ade a 
CUNY O rm ol ~ eee < ee 
Ree eal — on Ls GC: ao Zz cS — 
mal ko Y) {| ~~ a Cas — 4 ‘) 
Lo + “ee <r Bae) el et x 
ON eee a ee ava e ~f a _ x 
— YO amet Ot Ce a SS a > cr Oe 
ee = wend iY? ae ™N Belo eee é _ ~ 
i oa ome nw S3 ain ~ @) ~ “C> “C> is e — a 
co ee A ~ © a = oe) pO eee nd ok -_ pent ° -~ —~ 3% 
sm 2 oe << Co © < N WO ee) =~. > | ~ mG SS 
nom Le NY a =O) << = cS ae ~ re ers ae ft —<— mM) =f aS 
mn bee BOY ~ =a J a) ae a eS rt ee, pe ME | ca aes — . een "Sr 
lies eC Up ame em OR ~ <{ » NM hos Ae) Sih Geue)s <Co — 6 — iD 7” ~ @ any 
aN~e® ao Oost Os a, a rm O e — ®™N SON a eee tel a ~ aC ce Tt es 
eee OY BL ew | Ca a aa) ic) <a =~ WO -» ~ = ie LL oe | at mw “3° ww -L) oan yp 
Cet Sa tS) bee mime M eM of a Oo J >, Pagy  h n  e S  P fe es ees m= e a ~ ae 2 
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